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ABSTRACT 

Warm Dark Matter (WDM) has been invoked to resolve apparent conflicts of Cold 
Dark Matter (CDM) models with observations on subgalactic scales. In this work we 
provide a new and independent lower limit for the WDM particle mass (e.g. sterile 
neutrino) through the analysis of image fluxes in gravitationally lensed QSOs. 

Starting from a theoretical unperturbed cusp configuration we analyze the effects 
of intergalactic haloes in modifying the fluxes of QSO multiple images, giving rise 
to the so-called anomalous flux ratio. We found that the global effect of such haloes 
strongly depends on their mass/abundance ratio and it is maximized for haloes in the 
mass range 10 6 — 10 s M Q . 

This result opens up a new possibility to constrain CDM predictions on small 
scales and test different warm candidates, since free streaming of warm dark matter 
particles can considerably dampen the matter power spectrum in this mass range. As 
a consequence, while a (A) CDM model is able to produce flux anomalies at a level 
similar to those observed, a WDM model, with an insufficiently massive particle, fails 
to reproduce the observational evidences. 

Our analysis suggests a lower limit of a few keV (m v ~ 10) for the mass of warm 
dark matter candidates in the form of a sterile neutrino. This result makes sterile 
neutrino Warm Dark Matter less attractive as an alternative to Cold Dark Matter, in 
good agreement with previous findings from Lyman-a forest and Cosmic Microwave 
Background analysis. 
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1 INTRODUCTION 

The Cold Dark Matter (CDM) model has been successful 
in explaining a large variety of observational results such 
as the large scale structure of the universe and fluctua- 
tions of the Cosmic Microwave Background (CMB, Spergel 
et al. 2003, 2006). However, the CDM model faces some 
apparent problems on small scales: namely the overpredic- 
tion of galactic satellites, the cuspiness and high density of 
galactic cores and the large number of galaxies filling voids 
(Klypin et al. 1999, Moore et al. 1999a,b, Bode, Ostriker 
& Turok 2001, Avila-Reese et al. 2001, Peebles 2001 and 
references therein). These problems may well have complex 
astrophysical solutions. For instance the excess of galactic 
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satellites can be alleviated by feedback processes such as 
heating and supernova winds that can inhibit the star for- 
mation in low-mass haloes (Bullock, Kravtsov & Weinberg 
2001). 

Another natural cosmological solution to these prob- 
lems is to replace cold dark matter with a warm species 
(AWDM, see Bode, Ostriker & Turok 2001 and references 
therein). The warm component acts to reduce the small- 
scale power, resulting in fewer galactic subhaloes and lower 
central halo densities. 

One of the most promising WDM candidates is a sterile 
(right-handed) neutrino with a mass in the keV range; such 
a particle may occur naturally within extensions to the stan- 
dard model of particle physics (Dodelson & Widrow 1994, 
Dolgov & Hansen 2002, Asaka et al. 2005, Viel et al. 2005). 
A sterile neutrino is non-thermal in extensions of the mini- 
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mal standard model, with a life-time longer than the age of 
the universe. 

A strong constraint on the mass of WDM candidates 
comes from Lyman-a forest observations (neutral hydrogen 
absorption in the spectra of distant quasars), since they are 
a powerful tool for constraining the matter power spectrum 
over a large range of redshifts down to small scales. Recent 
analysis of SDSS quasar spectra combined with CMB and 
galaxy clustering data have set a lower limit on the mass 
of the sterile neutrino around m v ~ 10 — 13 keV (Seljak 
et al. 2006, Viel et al. 2006). In this paper we use a com- 
pletely different approach to put independent constraints 
on m v , using QSO gravitational lensing and the so-called 
anomalous flux ratio. 

Standard lens models, although they reproduce in 
general the relative positions of the images quite accu- 
rately, often have difficulties explaining the relative fluxes 
of multiply-imaged sources (Mao & Schneider 1998, Met- 
calf & Madau 2001, Dalai & Kochanek 2002, Metcalf and 
Zhao 2002), giving rise to the so-called anomalous flux ratio 
problem. 

Several possible explanations have been considered in 
the literature, the most plausible being that the lensing po- 
tential of real galaxies are not fully described by the sim- 
ple lens models used to compute lens characteristics. The 
most often invoked solution is to consider additional small- 
scale perturbations (i.e. dark matter haloes), which if lo- 
cated near a photon's light path can modify the overall lens 
potential (e.g. Raychaudhury et al. 2000, Saha et al. 2007) 
and significantly alter the observed flux ratio between dif- 
ferent images, in particular in the cusp or fold configuration 
(Metcalf & Madau 2001, Chiba 2002, Chen et al. 2003, Met- 
calf 2005a,b, Dobler & Keeton 2006). Those perturbers can 
be roughly divided in two categories: haloes that are inside 
the primary lens, usually referred as sub-haloes, and haloes 
that are along the line of sight, in between the source and 
the observer. This first category of haloes has been exten- 
sively studied in the past years both through analytic cal- 
culation (Metcalf & Madau 2001, Dalai & Kochanek 2002, 
Metcalf and Zhao 2002, Keeton 2003) and using numerical 
simulations (Bradac et al. 2002, Amara et al. 2006, Maccio 
et al. 2006). The latter two studies have came to the con- 
clusion that the impact of sub-haloes on lensing in the mass 
range 10 7 — 10 10 /i _1 Mq is very small. Even considering the 
impact of less massive subhaloes, usually not resolved in 
Nbody/hydro simulations, does not help in reproducing the 
observed number of anomalous flux ratios (Maccio & Mi- 
randa 2006). 

The effect of the second category of haloes, those along 
the line of sight, is still somewhat controversial (Chen 
et al. 2003, Metcalf 2005a,b). In particular Metcalf (2005a,b) 
found that dark matter haloes with masses around 10 6 — 
10 8 Mq can produce anomalies in the flux ratios at a level 
similar to those that are observed. The presence of a WDM 
particle even with a mass around 10 keV will strongly re- 
duce the number density of such small mass haloes, giving 
a different signature to the image fluxes. As a consequence, 
the observed anomalous flux ratios can be used to constrain 
the abundance of small haloes along the line of sight and 
therefore to put an independent constraint on the mass of 
the sterile neutrino as a possible WDM candidate. 

In this paper we analyze in detail the effect of subhaloes 



along the line of sight on an unperturbed cusp configuration 
in a ACDM model and in AWDM models with different val- 
ues of m v . We found that WDM models with m v < 10 keV 
fail to reproduce the observed anomalies in the lensed QSO 
flux ratios. Our results provide a new and independent con- 
straint on the mass of sterile neutrino, and they are in good 
agreement with previous constraints coming from Lyman-a 
forest and CMB analysis. 

The format of the paper is as follows: in section 2 we 
compute the expected halo abundance in different models; in 
section 3 we review briefly the lensing formalism we adopt. 
Section 4 is devoted to the description of our lensing simula- 
tions. In section 5 we present the numerical results, matching 
them with observations. We conclude with a short summary 
and discussion of our results in section 6. 



2 INTERGALACTIC HALO MASS FUNCTION 

The main goal of this work is to study the effect of dark mat- 
ter haloes along the line of sight on fluxes of QSO multiple 
images. In order to achieve it we first computed the number 
density of those haloes in the light cone between the source 
plane and the observer. 

For this purpose we used the Sheth and Tormen mass 
function (ST: Sheth & Tormen 2002), taking into account its 
evolution with redshift. We adopted a WMAPl-like cosmol- 
ogy (Spergel et al. 2003) with the following values for dark 
energy and dark matter density, normalization and slope of 
the matter power spectrum: Qa = 0.74, fi m = 0.26, <t$ = 0.9 
and n = l. 

The transfer function for the CDM model has been gen- 
erated using the public code CMBFAST (Seljak & Zaldar- 
riaga 1996). To compute the transfer function for WDM 
models we used the fitting formula suggested by Bode, Turok 
and Ostriker (2001): 

pWDM 

T2 ( k ) = ^CDM={ 1 + ( ak fT li,/ " (1) 

where a, the scale of the break, is a function of the WDM 
parameters, while the index v is fixed. Viel et al. (2005, see 
also Hansen et al. 2002), using a Boltzmann code simulation, 
found that v = 1.12 is the best fit for k < 5 h Mpc -1 , and 
they obtained the following expression for a: 

— (i&r (&r &r <*> 

This expression applies only to the case of thermal relics. 
In order to apply it to a sterile neutrino we take advan- 
tage of the one-to-one correspondence between the masses 
of thermal WDM particles (m x ) and sterile neutrinos (ra„) 
for which the effect on the matter distribution and thus 
the transfer function for both models are identical (Colombi 
et al. 1996). We used the m x — ra„ relation given by Viel 
et al. (2005), that reads: 

, ,„ /m^hc^aA 4 / 3 / 0.25^/3 /0.7\ 2 / 3 
m„ stcrilc = 4.43(^ 1 ^j ( — j (-) keV.(3) 

We used the expression given in eq:2 for the damping of 
the power-spectrum for simplicity and generality. More ac- 
curate expressions for the damping for concrete models of 
sterile neutrinos exist (Abazajian 2006, Asaka et al. 2007) 
and show that the damping depends on the detailed physics 
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Figure 1. Effects of WDM particles on the dark matter halo 
mass function at redshift zero. 



Figure 2. Evolution with redshift of the number of haloes above 
a fixed mass threshold in different models. The upper-most solid 
(blue) line is for M > 10 6 h~ 1 MQ in the ACDM model; the 
dashed and the dotted lines are for the same mass threshold but 
for a WDM mass of m„ = 10, 5keV respectively. The second set 
of (red) lines refers to a mass threshold of M > 10 7 /i -1 Mq. 



of the early universe in a rather non-trivial way. Naturally 
the results of this paper can be repeated using other expres- 
sions for the damping. 

The main effect of WDM is to dampen the power spec- 
trum of fluctuation on small scales, reducing the number of 
haloes at low masses (Bode, Turok & Ostriker 2001, Barkana 
et al. 2001, Paduroiu et al. 2007 in prep.). Figure 1 shows 
the ration between halo number density in WDM and CDM 
models as a function of the WDM mass m v . 

Typically lensed QSOs are located at a redshift around 
3. This implies that we also need to take into account the 
redshift evolution of the mass function in different models. 
Figure 2 shows the number of haloes more massive than 
10 6 /i _1 Mq (upper solid curve) and 10 7 ^"'Mq (lower solid 
curve) per Mpc cube at different redshifts. It is interesting 
to note that on such small mass scales the halo number den- 
sity tends to increase towards high redshift. We found that 
the evolution of the mass function, both in CDM and WDM 
models, can be well represented by the following fitting for- 
mula: 



logiV(> M,z) = iV +0.11 • 2° 



(4) 



where iVo is the logarithm of the halo number density at red- 
shift zero (No = log N(> M, z = 0)). The use of this fitting 
formula has the advantage of speeding up the calculation of 
the number of haloes in each lensing plane (see section 4). 

To conclude this section we want to emphasize that 
our particular choice of cosmological parameters does not 
influence the results we will present in the next section. 
For instance on the mass scales we are interested in (M < 
10 10 Ii~ 1 Mq) changing as from 0.9 to 0.7 would increase the 
number of haloes only by a few percent. 



3 LENSING FORMALISM 

We briefly recall the general expressions for gravitational 
lensing and refer, e.g., to the book by Schneider et al. (1992) 
for more details. The lens equation is defined as: 



9 = (3 + a(9) 



(5) 



where 0(9) is the source position and 9 the image position. 
d(0) is the deflection angle, which depends on k(6) the di- 
mensionless surface mass density (or convergence) in units 
of the critical surface mass density E cr i t , defined as: 



As 



(6) 



4ttG DlDls' 

where D$ , Dl, Dls ar e the angular diameter distances be- 
tween observer and source, observer and lens, source and 
lens, respectively. 

3.1 The cusp relation 

There are basically three configurations of four-image sys- 
tems: fold, cusp, and cross (Schneider & Weiss 1992). In 
this paper we will mainly concentrate on the cusp config- 
uration, that corresponds to a source located close to the 
cusp of the inner caustic curve (see figure 3). The behavior 
of gravitational lens mapping near a cusp was first studied 
by Blandford & Narayan (1986), Schneider & Weiss (1992), 
Mao (1992) and Zakharov (1995), who investigated the mag- 
nification properties of cusp images and concluded that the 
sum of the signed magnification factors of the three merg- 
ing images approaches zero as the source moves towards the 
cusp. In other words: 

fJiA + [i-B + /J-c 



Reus 



\ha\ + + Imc| 



0, for fitot 



(7) 
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Figure 3. Unperturbed cusp configuration: R c 



0.09. The 



source and image positions are marked by a solid circle and open 
triangles respectively. The opening angle is also shown. 



where ptot is the unsigned sum of magnifications of all four 
images, and A,B & C are the triplet of images forming 
the smallest opening angle (see figure 3). By opening an- 
gle, we mean the angle measured from the galaxy center 
and spanned by two images of equal parity. The third im- 
age lies inside such an angle. This is an asymptotic relation 
and holds when the source approaches the cusp from inside 
the inner caustic " astroid" . This can be shown by expand- 
ing the lensing map to third order in the angular separation 
from a cusp (Schneider & Weiss 1992). Structure on scales 
smaller than the image separation will cause R cusp to differ 
from zero fairly independently of the form of the rest of the 
lens. Note that by definition of R CU s P used here, it can be 
either positive or negative. A perturber is more likely to re- 
duce the absolute magnification for negative magnification 
images (Metcalf & Madau 2001, Schechter & Wambsganss 
2002, Keeton et al. 2003) and to increase it for positive par- 
ity images. As a result, the probability distribution of R CU s P 
will be skewed toward positive values. 



3.2 The unperturbed lens 

We used the GRAVLENS code (Keeton 2001)* to create 
a lens configuration for which the cusp relation is roughly 
satisfied (see figure 3). The main, smooth, lens has been 
modelled as a singular isothermal ellipsoid (SIE) (Kormann, 
Schneider, & Bartelmann 1994) to take advantage of its sim- 
plicity. This model has been widely used in lens modeling 
and successfully reproduces many lens systems (e.g. Keeton 
et al. 1998, Chiba 2002, Treu & Koopmans 2002). The el- 
lipsoidal primary lens has a mass equal to 5 x 10 11 Mq, is 
oriented with the major axis along the y axis in the lens 



t The software is available via the web site http://cfa- 
www.harvard.edu/castles 



plane and has an ellipticity of 0.33. The redshift of the lens 
has been fixed to zi = 0.3 in agreement with typical ob- 
served ones (i.e. Tonry 1998). The cusp relation, defined by 
equation 7, for this smooth lens gives R CU s P = 0.09, and this 
is one of the configurations previously studied in Maccio and 
Miranda (2006, namely Config2). We tested that our results 
do not depend on this particular choice for the unperturbed 
configuration and do apply to any cusp configuration. 



4 SUBHALOES ALONG THE LINE OF SIGHT: 
IDEA AND PROCEDURE 

The purpose of this work is to compute the effects of inter- 
galactic haloes, along the line of sight, on an unperturbed 
cusp lensing configuration to extract information on the 
matter power spectrum on small scales. In this approach, 
we model our haloes as singular isothermal spheres (SIS). 
A SIS, with density profile p tx r~ 2 , is a simple model that 
is often used in lensing because its simplicity permits de- 
tailed analytic treatment (e.g., Finch et al. 2002). The model 
has been used to represent mass clumps for studies of sub- 
structure lensing, after taking into account tidal stripping by 
the parent halo (Metcalf & Madau 2001; Dalai & Kochanek 
2002). Again, the simplicity of the SIS makes it attractive for 
theoretical studies: a tool that not only reveals, but also elu- 
cidates, some interesting general principles. For the 1O 6 M0 
haloes relevant for this work, the SIS profile does not differ 
dramatically from the NFW (Navarro, Frenk, & White 1996) 
profile inferred from cosmological N-body simulations (Kee- 
ton 2003). Besides, the SIS model yields conservative results. 
Since an NFW halo is centrally less centrally concentrated 
than a SIS halo, it is less efficient as a lens and therefore 
would have to be more massive in order to produce a given 
magnification perturbation. Maccio & Miranda (2006) have 
shown that a SIS model will induce lensing effects marginally 
stronger then those caused by an NFW profile with concen- 
tration parameter c ~ 55, corresponding to a mass around 
10 6 Mq. Haloes in a WDM model are expected to be less 
concentrated due to the top-down structure formation sce- 
nario (Eke, Navarro & Steinmetz 2001, Paduroiu et al. 2007 
in prep). In this case the SIS approximation can possibly 
overestimate the total effect of WDM perturbers, making 
our lower bound to the WDM particle mass even stronger. 

A SIS halo model is completely characterized by its Ein- 
stein radius: 



0e 



47rq 2 D LS 
c 2 D s '' 



(8) 



where a is the halo velocity dispersion, and Ds,Dls are the 
angular diameter distances introduced in sec. 3. We adopt a 
source redshift z s — 2. 

We filled the portion of Universe along the line of sight 
with cubes, then the subhaloes inside each cube were pro- 
jected onto the middle plane (see figure 4) . We used a total 
of 100 different lens planes roughly equally distributed in 
space between the source and the observer. This results in 
iVi = 85 planes behind the main lens and N2 = 15 planes in 
front of it. The size of the cubes was defined as follows. 

Two close planes were separated by Azi = (z max — 
Zi)/Ni if situated behind the main lensing galaxy, and by 
AZ2 = (zi — z m i n )/N2 for planes in front if it, where z m in = 
0.01 and z max = z s — 0.1. 
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Figure 4. A schematic diagram of the type of lensing system 
being considered. There is one primary lens responsible for the 
multiple images of the source. In addition, there are many sec- 
ondary lenses (most not shown). The unperturbed light paths are 
deflected only by the primary lens and with an appropriate model 
for the primary lens will meet on the source plane. If the deflec- 
tions from secondary lens planes arc taken into account without 
changing the primary lens model, the light will follow the per- 
turbed light paths (dashed curves). This diagram is not to scale 
in any respect. 



The size of a comoving volume inside a solid angle dQ, 
and a redshift interval dz is given by (Hogg 1999): 



dV c = Dh 



E(z) 



dQ dz 



(9) 



where Da is the angular diameter distance at redshift z and 
E(z) is defined as: 



E(z) = yjQ M (1 + «) 3 + fifc (1 + z) 2 + Oa 



(10) 



with Qm, f2 fe and FIa the density parameters of matter (cold 
and warm), curvature, and cosmological constant, respec- 
tively. 

We populated each cube with dark matter haloes, whose 
total number and mass distribution was chosen according to 
the ST mass function at the appropriate redshift (see sec- 
tion 2). Halo positions and redshifts (within Azi^) were ran- 
domly assigned. Within a solid angle dfi of 3" x 3" squared 
arcsec, the total number of haloes with mass larger than 
10 6 M Q comes to 512 for the ACDM model adopted in this 
paper. This number drops in a consistent way in a warm 
dark matter scenario, depending on m„. For a WDM parti- 
cles mass of 10 keV we obtain 238 haloes along the line of 
sight within the same dfl, and even fewer (156, 135) for a 
less massive choice for m v (7.5, 5 keV, see fig 1). 

Since we are interested in flux anomalies, we consider 
only cases in which we do not have image splitting due to the 
extra haloes along the line of sight. Therefore we do not allow 
any of those haloes to be closer than twice its Einstein radius 
(9e) from any images, in order to prevent image splitting 
(see Schneider, Ehlers & Falco 1992 and references therein). 
On average only few haloes ( 3, for LCDM) fail in satisfy- 
ing this criterium and we tested their removal/inclusion do 
not affect the final R distribution in any way. Let r\ denote 
the two-dimensional position of the unperturbed image with 
respect to the perturber on the I plane , measured with re- 
spect to the intersection point of the optical axis with the / 
plane and £ the light ray impact parameter on the 7' plane. 



In the absence of image splitting a SIS perturber will affect 
the position of each image according the following: 



r 1 = tjp--a{t)D 1 , I . 



(11) 



Introducing the angular coordinates r\ — DjQi and £ = 
D'jO'j, and given that a(£) = 9e for a SIS, the equation 
for the flux becomes 



A* = 



(12) 



where the quantities with subindex I refer to the (unper- 
turbed) image position with respect to the perturber and so 
Di, T>ii , Djij are the distances between observer and the / 
plane, observer and I' plane, / plane and I' plane, respec- 
tively. On each single lens plane the total effect on the image 
magnification factor /i is obtained by summing up contribu- 
tions by each perturber. In principle one should sum the 
magnification tensors first and then take the determinant. 
The two methods (scalar or matrix sum) do not lead to the 
same result because det(A + B) ^ det(A) + det(B). In the 
case of scalar sum and two SIS perturbers with Einstein radii 
8e,i and 9e, 2, the total magnification depends on the order 
in which the two lenses act on the source: /ii^ is different 
from H2,i- The error introduced by a direct sum is of the 
order of the ratio between the /xi ; 2 and ^12,1 • This quantity 
can be directly computed from eqs: 11 and 12 and it is al- 
ways < max(#E,i, 0e,2)//3. In our case, due to the low mass 
of our perturbers, the ratio 0i/9e,i is of the order of 200- 
800, which gives an error less than 1% for the total fi. There 
is still a small chance to have a substructure located at a 
place where 61 ~ 6e,z- We looked for this possibility and it 
happened only 8 times over 100.000 substructure position 
realizations, giving a negligible effect on the final averaged 
value of Rcusp. 

Generally a matter clump will change the positions of 
the images slightly, so if a lens model is chosen to fit the ob- 
served image positions perfectly it will not do it anymore af- 
ter the perturber is added. To produce a perfectly consistent 
lens model one would have to adjust the main lens model 
for each realization of the intergalactic haloes. This is very 
computationally expensive and not necessary in practice. 
The shifts in positions are generally small when the masses 
of the secondary lenses are small (f» 0.1" for M ~ 10 s Mq 
Metcalf 2005a) and, in addition, since the host lens model 
is degenerate it is ambiguous how it should be adjusted to 
correct for the shift. The goal here is to reproduce all the sig- 
nificant characteristics of the effects induced by the observed 
lens (image configuration, fluxes) so that one can determine 
whether lenses, that look like the observed ones and have the 
observed ratio anomalies, are common in CDM/WDM mod- 
els. For the source, we adopt the point-like approximation. 
The importance of considering the source size lies mainly in 
the capability to disentangle different subhaloes mass limits 
(Chiba et al. 2005, Dobler & Keeton2006). As remarked by 
Chang & Refsdal (1979) and many authors afterwords (see 
Metcalf 2004 and references therein), the projected size (on 
the lens plane) of the emitting regions of QSOs are expected 
to be different and this can be used to remove, eventually, 
lens model degeneracy and improve the sensitivity to sub- 
structure properties. In our cases, the size of the radio emit- 
ting region, when projected on the lens plane, is expected 
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to be affected by structures with masses larger than 10 5 Mq 
(Metcalf 2005a,b). 

In a single realization of our perturbed lens configura- 
tion the light coming from the source is deflected by w 500 
haloes (plus the main lens) before reaching the observer. 
Each one of the three images forming the cusp configu- 
ration is shifted and amplified, giving as a result a modi- 
fied Rcusp value, different from the original (unperturbed) 
one of Rcusp = 0.09. Sometimes, when a massive halo 
(M > 10 8 Mq) happens to be close to one of the images, 
this image can be strongly deflected, resulting in a breaking 
of the cusp configuration. In the statistical studies presented 
here these cases are simply excluded from the final sample. 
In total we performed 2,000 realizations (with different ran- 
dom seeds for generating masses and positions of perturbers) 
of each model (CDM/WDM), obtaining 2,000 different final 
lensing configurations. For some of these final configurations 
(with high Rcusp values), we try to fit image positions and 
magnification factors with the GRAVLENS code, using a 
smooth lens model. While is relatively simple to reproduce 
the image geometrical properties, it is never possible to get 
the right flux ratios, with such a simple model. 



5 RESULTS 

The first part of this section is devoted to presenting the 
effects of haloes along the line of sight (l.o.s.) on the cusp 
relation in a standard (A)CDM scenario. The plots show the 
probability distribution for the cusp relation value, consid- 
ering 2,000 different realizations of the same model. Those 
realizations share the same total number of perturbers, but 
differ in their masses (randomly drawn from a ST distribu- 
tion), positions (randomly assigned within the lens plane) 
and redshifts (randomly chosen within Azi^)- 

The cusp relation defined by equation 7 holds when the 
source is close to the cusp. As soon as the source moves away 
from the cusp, deviations from Rcusp = are observed, even 
for the smooth lens model. On the other hand the closer the 
source is to the cusp, the smaller is the angle spanned from 
the three images. Therefore, in order to take into account 
the position of the source in evaluating the cusp relation, it 
is better to define the anomalous flux ratio as: 

R = Rcusp (13) 

where A6 is the opening angle spanned by the two images 
with positive parity defined from the center of the galaxy. 
With this new definition of the cusp relation a set of three 
images is said to violate the cusp relation if R > 1. This 
makes the comparison between simulations and observations 
much more straightforward. For this comparison we used 
the same data presented in Maccio et al. (2005) . There are 
5 observed cusp caustic lenses systems ( summarized in ta- 
ble 1): B0712+472 (Jackson et al. 1998), B2045+265 (Koop- 
mans et al. 2003), B1422+231 (Patnaik & Narasimha 2001), 
RXJ1131-1231 (Sluse et al. 2003) and RXJ0911+0551 (Kee- 
ton et al. 2003); the first three are observed in the radio 
band, the last two in optical and IR. Three of them violate 
the reduced cusp relation (i.e. R > 2n/A9). 

Figure 5 shows the 7? probability distribution for the 
three possible categories of perturbers. The dotted (red) 
line shows the effect of subhaloes inside the primary lens 



lens 


A6> 


7?cusp 


obs. band 


B0712+472 


79.8° 


0.26 ± 0.02 


radio 


B2045+265 


35.3° 


0.501 ± 0.035 


radio 


B1422+231 


74.9° 


0.187 ±0.006 


radio 


RXJ1131-1231 


69.0° 


0.355 ± 0.015 


optical/IR 


RXJ0911+0551 


69.6° 


0.192 ± 0.011 


optical/IR 



Table 1. The image opening angles and cusp caustic parameters 
for the observed cusp caustic lenses. 



that can be directly tested by current numerical simulations 
(i.e with masses > 10 7 M , Maccio et al. 2006). The short- 
dashed (cyan) line shows the effect of lower mass subclumps 
(still inside the primary lens) as measured by Maccio and 
Miranda (2006). The solid (blue) line shows the effect of the 
haloes along the line of sight considered in this work; here 
we considered only haloes with M > 5 x 10 6 Mq . As already 
noticed the first two categories of perturbers fail in repro- 
ducing the high value tail that arises in the observational 
data around R = 2. On the contrary, the signal coming from 
haloes along the l.o.s. has a probability distribution which 
remains almost flat in R range 1-2, where 2 (out of 5) of the 
observed systems lay. 

Thanks to this pronounced tail at high R value, haloes 
filling the light cone between the source and the observer can 
easily account for all the observed cusp systems, providing 
a solution to the anomalous flux ratio issue. Our results are 
in fair agreement with those previously obtained by Metcalf 
(2005b) and seem to confirm that a previous result on the 
same subject obtained by Chen et al. (2003) did underesti- 
mate the effects of intergalactic structure. Chen et al. (2003) 
used the cross section (or optical-depth) method to calcu- 
late the magnification probability distribution. This method 
is mainly valid for rare events and this is not the case since, 
as shown in section 4, the number of lensing events is of 
the order of 500. A more detailed and general comparison of 
the two methods can be found in Metcalf (2005b). In Met- 
calf (2005b) the author used an approach similar to ours 
making a direct lensing simulation in order to compute the 
effects of haloes along the l.o.s., modelling them using an 
NFW density profile. Although in his work the author ana- 
lyzed each observed configuration separately, finding slightly 
different individual R probabilities for different systems, the 
similarity of the results is a good proof a posteriori that 
our assumptions of SIS parametrization for perturbers and 
point-source approximation did not introduce a strong bias 
in the results. 

In the previous analysis we restricted the mass range 
to haloes more massive than M — 5 x 10 6 Mq. In figure 6 
the probability distribution for R is shown for two differ- 
ent choices of the minimum halo mass: M > 5 x 10 6 Mq 
(solid, blue line) and M > 10 5 Mq. In the latter case the 
total number of structures is around 5,500 and the lensing 
simulation code slows down considerably. A close compar- 
ison of the two histograms clearly shows that considering 
less massive haloes does not not improve the results sub- 
stantially; so in the following we will only consider haloes 
with M > 5 x 10 6 M . 

In some cases, when the averaging process is restricted 
to a lower number of realizations (~ 200) we found that the 
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Figure 5. R probability distribution for different categories of 
(sub)haloes within the CDM scenario. The dotted line shows the 
effect of substructures (with M > 10 7 Mg) inside the lens galaxy 
(Maccio et al. 2006); the long-dashed line is for less massive sub- 
haloes (M = 10 5 - 10 7 M ) still inside the primary lens (Maccio 
& Miranda 2006). The solid line is for the haloes along the line of 
sight with mass > 5 X 10 6 Mq studied in this work. Observational 
results are also shown (long dash histogram). 



observational data are reproduced with a high confidence 
level as shown in figure 7. These results are probably due 
to effects induced by single massive perturbers close to a 
particular image: or a positive image is highly magnified or 
a negative one is demagnified (note that in equation 7 we 
consider the absolute values for fj,i), providing an anoma- 
lous R. While with a low number of realizations (~ 200) 
these single events contribute significantly to the global R, 
a higher number of realizations (> 10,000) permits all the 
images to be affected by massive clumps, smoothing the final 
probability distribution. 

The introduction of a WDM particle damps the mat- 
ter power spectrum on small scales, reducing the number 
of haloes along the l.o.s. In figure 8 we show the probability 
distribution of R as a function of the mass of the WDM can- 
didate. Changing the WDM particle mass from ra„ = 12.5 
to m v — 7.5 keV drops the tail at R — 2 from a 10% proba- 
bility to a 1.5% one. For m„ = 5 keV we have a P(R) higher 
than 5% only for R < 1.3. In the latter case only 20 haloes 
are inside the volume sampled by the three images, and this 
model tends to leave the value of R close to the unperturbed 
one. A model with a 10 keV sterile neutrino, if compared to 
a model with m v = 12.5 keV, gives a slightly lower proba- 
bility (8% vs 10%) to have a configuration with R — 2. Due 
to the limited number of observed cusp systems it is hard to 
disentangle those two models, and we think that it is fair to 
say that m„ = 10 keV is still in agreement with the data. 

Figure 9 shows the comparison between the observa- 
tional data, the standard (A) CDM model and a WDM model 
with a sterile neutrino mass of 12.5 keV, which is close to the 
current limit provided by Lyman-a + CMB analysis (Seljak 



Figure 6. R distribution for haloes along the line of sight for 
two choices of their minimum mass: M > 10 5 (dot line) and 
M > 5 X 10 6 Mq (solid line). The dashed histogram shows the 
observational data. 




Figure 7. R probability distribution for CDM considering a lower 
number of realization (ps 200) in the averaging process (see text). 
The dashed histogram shows the observational data. 



et al. 2006). In this case in both the warm and cold dark 
matter scenario, haloes along the line of sight can easily ac- 
count for the two observed cusp systems with R ss 2, offering 
a viable solution to the anomalous flux ratio issue. On the 
contrary a warm dark matter model with less massive par- 
ticles (i.e. with a higher free streaming scale length) fails 
in reproducing the observational data due to the reduced 
number density of haloes along the line of sight. 
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Figure 8. Probability distributions for different warm particle 
masses: m„ = 5keV (long-dashed line), m„ = 7.5 keV (dot line), 
m u = 10 keV (dashed-dot line), m u = 12.5 keV (solid line). 
Dashed line shows the probability distribution of observational 
data. 



I " " I " " I " " I " " I " " I " " I " " I " " I 
WDM vs CDM, 5xl0 6 -10 9 M 



m„ = 12.5 keV 




Figure 9. R distribution probability for: observed values (dashed 
line), CDM haloes more massive than 5 X 10 6 A/q (solid line) and 
WDM subhaloes with m v =12.5keV (dotted line). 



6 DISCUSSION AND CONCLUSIONS 

Interest in warm dark matter models has been sporadic over 
the years, although this class of models could help allevi- 
ate several problems on small scales that occur with cold 
dark matter. In order to constrain the WDM scenario, pre- 
cise measurements of the matter power spectrum on small 
scales are needed; for this purpose Lyman-a forest and CMB 



data have been extensively used (Seljak et al. 2006, Viel 
et al. 2006). 

In this paper we show that image flux ratios in multiple 
gravitationally lensed QSOs can be modified by haloes along 
the line of sight in the mass range 10 6 — 10 7 Mq; this effect 
opens a new window to study the matter power spectrum 
on small scales and provides a new and independent method 
to constrain the mass of WDM candidates m v . 

The observed anomalous flux ratio in lensed QSOs can 
be explained by adding small perturbations to the smooth 
model used to parametrize the main lenses. Those per- 
turbers can be identified with dark matter haloes that 
happen to be close to the images' light paths. Recent 
results based on numerical N-Body (Amara et al. 2006, 
Rozo et al. 2006) and hydrodynamical simulations (Maccio 
et al. 2006) have shown that it is hard to reconcile the ob- 
served high number of cusp relation violations with the total 
number of substructures inside the primary lens predicted 
by the ACDM model. This is true even when the limited 
mass resolution of numerical simulations is taken into ac- 
count (Maccio and Miranda 2006). 

The hierarchical formation scenario predicts that the 
universe should be filled by a large number (more then 
10 3 per fo -1 Mpc 3 ) of dark matter haloes with masses 
M « 10 6 M . We employed the Sheth & Tormen mass func- 
tion to estimate the expected number of haloes in this mass 
range along the line of sight of lensed QSOs. We found that 
on average there are more than 500 haloes in between the 
source and the observer, within a light cone with an aperture 
of 3 arcsec. Using direct lensing simulations and a singular 
isothermal sphere approximation we computed the effects of 
those haloes on an unperturbed cusp configuration. We gen- 
erated more than 10 4 different realizations of our global (lens 
+ perturbers) lensing system, varying masses, positions, and 
number of haloes. 

We found that on a statistical basis (averaging on dif- 
ferent realizations) this class of perturbers can modify con- 
sistently the fluxes of QSO multiple images at a level com- 
parable to the observed one, in good agreement with previ- 
ous studies on this subject (Metcalf 2005a,b). In some cases 
when the averaging process is restricted to a lower number 
of realizations (« 200, see figure 7) we found that the obser- 
vational data are reproduced with a high confidence level. 

An important result of our study is that the bulk of 
the signal on QSO fluxes is due to haloes in the mass range 
10 6 - 10 7 M Q . Since the number density of such haloes, and 
therefore their effect on the cusp relation, can be strongly 
damped by the presence of a WDM candidate, the observed 
number of anomalous flux ratios can be used to constrain 
the mass of WDM particles. 

Adding an exponential cut-off to the transfer function 
of WDM models we computed the number density of small 
haloes as a function of the mass of the warm particles. We 
show that if WDM is due to a sterile neutrino, then, in 
models with m v < 10 keV, the number of dark haloes along 
the line of sight is too low to affect in a consistent way 
the fluxes of lensed QSOs, failing to reproduce the observed 
abundance of systems with high R values. This lower limit 
for the mass of the sterile neutrino is in good agreement with 
results obtained using different methods. 

The main limitation of this study is represented by the 
few observational data that are available in the literature. 
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However, future experiments such as Dune, are likely to ob- 
serve more then 1000 lensed quasars, of which several hun- 
dreds should be quadruples due to the magnification bias. 
It will provide new lensing systems to be analyzed and thus 
more tightly constrain the WDM scenario. 
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